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Abstract. We present a theoretical model for the non-thermal radio emission from single hot stars, in terms of synchrotron 
radiation from electrons accelerated in wind-embedded shocks. The model is described by five independent parameters each 
with a straightforward physical interpretation. Applying the model to a high-quality observation of Cyg OB2 No. 9 (05 If), 
we obtain meaningful constraints on most parameters. The most important result is that the outer boundary of the synchrotron 
emission region must lie between 500 and 2200 stellar radii. This means that shocks must persist up to that distance. We also find 
that relatively weak shocks (with a compression ratio < 3) are needed to produce the observed radio spectrum. These results 
are compatible with current hydrodynamical predictions. Most of our models also show a relativistic electron fraction that 
increases outwards. This points to an increasing efficiency of the acceleration mechanism, perhaps due to multiple acceleration, 
or an increase in the strength of the shocks. Implications of our results for non-thermal X-ray emission are discussed. 
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1. Introduction 

Many hot stars of spectral type O and B are observable at radio 
wavelengths, due to thermal emission from the circumstellar 
ionized gas in the stellar wind. The thermal radiation is free- 
free emission (Bremsstrahlung) from an electron accelerated in 
the Coulomb field of an ion. The emergent radio spectrum has a 
characteristic spectral index 1 a a 0.6 (Wright & Barlow 1975), 
which fits the observations for most O stars quite well. 

However, 25% of the brightest O stars have a radio spec- 
trum with a spectral index which is very different from the 
thermal wind emission (Abbott et al. 1984). It is believed 
that this non-thermal radio emission is synchrotron radiation 
by relativistic electrons (White 1985). The preferred mech- 
anism to accelerate electrons to relativistic energies is parti- 
cle acceleration in shocks (Fermi 1949), as shocks are known 
to exist in stellar winds. In single stars, shocks are gener- 
ated by the instability of the driving mechanism of the wind 
(Owocki & Rybicki 1984). For binaries, shocks also arise in 
colliding winds (Eichler & Usov 1993). Non-thermal radio 
emission from colliding wind binaries has been further inves- 
tigated by Dougherty et al. (2003). In this paper we limit our- 
selves to single stars. 

Besides shocks, a magnetic field is needed to produce a 
synchrotron spectrum. Magnetic fields of normal hot stars are 
below current detection limits (~ 100 Gauss). It is plausible 
that a small magnetic field is present in the star, either the fos- 
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sil magnetic field, or a field created by a dynamo mechanism 
(MacGregor & Cassinelli 2003). 

The model proposed by White (1985) was further devel- 
oped by Chen & White (1991, 1994), who showed that a syn- 
chrotron model can successfully reproduce the radio spectrum 
of 9 Sgr and Cyg OB2 No. 9. 

In the present paper, we investigate what properties a syn- 
chrotron model should have in order to reproduce an observed 
non-thermal radio spectrum. More specifically, we look at all 
models that fit a given set of observations and thus derive the 
permitted range of a number of parameters. The stringency of 
the constraint on the parameters depends critically on the qual- 
ity of the available observations, as well as on the number of 
observed wavelengths. We selected Cyg OB2 No. 9 as the ob- 
ject of this study because a set of three simultaneous observa- 
tions (at 2, 6 and 20 cm) with small error bars is available. 

For the sake of clarity it is desirable that the properties of 
the model be contained in a limited number of parameters with 
a clear physical interpretation. Rather than calculating the mo- 
mentum distribution of relativistic electrons from a model of 
the shocks in the stellar wind, we use a power-law dependence 
on the compression ratio. The uncertainty about the exact shock 
structure is so large that an ab initio calculation of the mo- 
mentum distribution would lead to a greater complexity of the 
model without increasing its physical accuracy. 

The remainder of this paper is organised as follows. In 
Sect. 2 we describe the model for the synchrotron radiation 
in detail. The effect of cooling mechanisms on the momentum 
distribution is investigated in Sect. 3. In Sect. 4 we apply the 
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model to Cyg OB2 No. 9. Finally, we discuss the results in 
Sect. 5 and give some conclusions. 

2. The model 

2.1. Power 

In the presence of shocks and a magnetic field, electrons are ac- 
celerated to high velocities (Bell 1978) through the first-order 
Fermi mechanism (Fermi 1949) and emit synchrotron radiation 
as they gyrate around the magnetic field lines at relativistic ve- 
locities. The synchrotron power radiated per unit volume per 
unit frequency by a single electron with momentum p, moving 
in a magnetic field B (dependent on distance r) with pitch angle 
6, can be expressed as (Westfold 1959) 



P v (p,6,r) 
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m e c 2 
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where m e and e are the electron mass and charge, c the speed 
of light, v the frequency, F(x) = x f dr/K^/j, (rf) and 
the modified Bessel function. For the meaning of f(v, p), see 
Sect. 2.2. The synchrotron power spectrum extends up to fre- 
quencies of order v s before falling away, with the critical fre- 
quency v s defined as 
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2.2. Razin effect 

The presence of f(v, p) in the synchrotron power is a conse- 
quence of the Tsytovitch-Razin effect (or Razin effect for short; 
Tsytovitch 1951; Razin 1960). Relativistic beaming plays a 
dominant role in the explanation of synchrotron emission. In 
the presence of a cold plasma, the refraction index is smaller 
than unity. This increases the beaming angle, thereby reduc- 
ing the beaming effect. Thus the synchrotron emitting power 
is greatly reduced. Ginzburg & Syrovatskii (1965) expressed 

f(y,p) as 



f(y,p) 
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with V() = y 2s£_ the plasma frequency and n e the number den- 
sity of thermal electrons. 

In the absence of the Razin effect (/(v, p) = 1), the essential 
part of the synchrotron emission of an individual electron with 
momentum p is emitted below the critical frequency v s , with a 
maximum at 



v = O.29v s sin0« 1.3 



Bsin9 
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MHz. 
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Away from the maximum, the synchrotron radiation falls off 
rapidly. We can therefore use Eq. (4) to estimate the momen- 
tum range relevant for the radio spectrum. For a magnetic field 



of 0. 1 Gauss, radiation between 2 - 20 cm ( 1 5 - 1 .4 GHz) comes 
from electrons with momenta between « 50 - 150 MeV/c, 
where we set sin 6 = 1 . 

From Eq. (3) we see that the Razin effect is important 
when v «: vop/(m e c). If on the other hand v » vop/(m e c), 
the influence of the Razin effect is small. We can then use 
Eq. (2) to substitute p and find the frequencies for which the 
influence of the Razin effect is small. In the literature (e.g. 
Ginzburg & Syrovatskii 1965), this is given by 

tie 

v » 20-^ (5) 

where, in view of the approximate nature, the factor 0.29 has 
not been taken into account. In a typical hot-star wind, the ther- 
mal electron density at 100 R„ is of order 10 7 cm 4 . For a mag- 
netic field of 0.1 Gauss this means that the Razin effect begins 
to suppress radiation at wavelengths larger than 15 cm, well 
within the observable range. For smaller magnetic fields, the 
radiation is suppressed at even smaller wavelengths e.g. 1.5 cm 
for a 0.01 Gauss field. 



2.3. Flux 

The synchrotron emissivity of a distribution N(p, r) of relativis- 
tic electrons can be expressed as 



1 r Pc 

Mr) = — dpN(p, r)P Y (p, r), 

4n Jpo 



(6) 



where P v (p, r) is the synchrotron power P v (p, 6, r) averaged 
over solid angle assuming that the electron velocity distribu- 
tion is locally isotropic. The low (high) momentum cut-off of 
the electron distribution is defined by po (respectively p c ). We 
will specify the cut-off values in Sect. 3. 
The emergent flux is given by 



dr4nr 2 j v (r), 



(7) 



where D is the distance to the star. The integration boundaries 
reflect the fact that, at any given radio frequency, only a lim- 
ited part of the wind contributes to the flux. The choice of these 
boundaries is rather subtle. The lower boundary mimics the ef- 
fect of the large free-free optical depth of the wind, which effec- 
tively shields the inner part of the synchrotron emitting region 
from observation. The lower boundary must increase with in- 
creasing wavelength, to reflect the increasing free-free optical 
depth of the wind. It is customary (e.g. White 1985) to use the 
characteristic radius from Wright & Barlow (1975) as the lower 
boundary. However, this radius corresponds to a small optical 
depth (t v « 0.25) and was introduced to produce the correct 
thermal radio flux with a simple one-dimensional integral. Its 
physical meaning is often misstated and is actually quite lim- 
ited. It is not true that all observed radio emission originates 
from above R v . It is quite easy to replace Eq. (7) by the correct 
expression which is, however, rather cumbersome. We there- 
fore relegate its derivation to Appendix A and prefer to use the 
simple approximation Eq. (7) in the discussion of our results. 
The calculations were done using the exact expression given in 
Eq. (A.6). 
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As for the outer integration boundary, it is well known 
that instability-generated shocks decay as they move out with 
the wind (e.g. Runacres & Owocki 2002). It can therefore be 
expected that they eventually become too weak to accelerate 
electrons to relativistic energies. Because relativistic electrons 
rapidly lose momentum as they move away from the shock 
(Sect. 3.4), the outer boundary of the synchrotron emission re- 
gion (R max ) can also be visualised as the position of the last 
shock strong enough to accelerate electrons to relativistic en- 
ergies. It does not depend on frequency and is a fundamental 
parameter of our models. 

The total emergent flux has a minor contribution from 
free-free emission by thermal electrons. This contribu- 
tion is calculated using the Wright & Barlow formalism 
(Wright & Barlow 1975) and added to the non-thermal flux. 

2.4. Number distribution and magnetic field 

In order to derive an explicit expression for the emissivity, we 
need to specify the number distribution N(p, r) and the mag- 
netic field B(r). 

It was shown by Bell (1978) that the acceleration of elec- 
trons through the Fermi mechanism results in a power-law mo- 
mentum distribution with exponent n, where n is a function of 
the compression ration of the shock: n = (x + 2) /Of - 1). For 
reasons of simplicity, we also assume that the distribution is 
a continuous function of distance, given by a power law with 
exponent 6. N(p, r) is then expressed as 



(8) 



N(p,r) = N \ — 



where is a normalisation constant related to the number of 
relativistic electrons. The total number of relativistic electrons 
at the stellar surface is given by f Pc dpNo p~ n . This gives 



No = f*n* e (n - \)pl 



n-l 



(9) 



where we introduced /„ and n* as the fraction of relativistic 
electrons and the total electron number density at the stellar 
surface. We also assumed p c » p Q (see Sect. 3). The radial 
variation of the fraction of relativistic electrons in the wind 
is described by f*(r / R t y s+2 . It must be emphasised that /» is 
merely a parameter used to set the overall level of this fraction. 
It should not be too strictly interpreted as a physical quantity. 
Indeed, close to the stellar surface, the relativistic electron pop- 
ulation may be strongly suppressed by inverse Compton cool- 
ing in the intense UV-radiation. Therefore f, need not represent 
the actual relativistic electron fraction at the stellar surface. 
For n* we use 




model parameters: 

n=2.85;5=1;R =716 R,; 

max 

B.=100 Gauss; f.=1. 02x10" 
p =1 MeV/c;p=15GeV/c 



I (cm) 

Fig. 1. The full line represents a model that fits the VLA-observations 
(21 Dec 1984) for Cyg OB2 No. 9 (Bieging et al. 1989), while the 
dashed line neglects the Razin effect. The diamonds are the observa- 
tions and the error bar is given by the size of the symbol. The stellar 
parameters are listed in Table 1. 



the velocity has not been taken into account), but merely pro- 
vides the correct scaling factor for the electron density at large 
distances. 

In Eq. (1) the magnetic field B must be specified. 
Unfortunately, the magnetic fields of hot stars are below cur- 
rent detection limits. Therefore, we can say little about their 
radial dependence. For want of anything better, we will adopt 
the expression (Weber & Davis 1967) 



Voo r 



(11) 



where B t is the surface magnetic field and v rot the rotational 
velocity of the star. This expression is only valid at large dis- 
tances from the star (r > 10 R*), which covers the region where 
the synchrotron radiation is emitted. For the surface mag- 
netic field, the detection limit, which is of order ~ 100 Gauss 
(Mathys 1999), serves as an upper limit. A tentative lower limit 
can be derived from the synchrotron emission itself. At local 
magnetic fields smaller than ~ 0.005 Gauss at 100 R„, all ra- 
dio emission would be suppressed by the Razin effect. Using 
Eq. (11) and v ro t/Vco = 250/2900, this translates to a surface 
magnetic field of order 5 Gauss. 

Once a choice for the parameters B„, n, 6, R mdx , po and 
p c has been made, Eq. (A. 6) can be integrated numerically. The 
total emergent flux (i.e. non-thermal + thermal) for one such 
model is shown in Fig. 1. 



n„ = 



yM 



AnRlvccpmn 



(10) 



where M is the mass loss rate, Voo the terminal velocity of the 
stellar wind, m H the mass of a hydrogen atom, y the mean num- 
ber of electrons per ion and p the mean ionic weight. For sim- 
plicity we assume y = 1 and p = 1 .4. Note that n* is not the 
actual electron density at the stellar surface (as the variation of is found for the synchrotron flux (Ginzburg & Ozernoy 1966), 



2.5. Tests 

If we neglect the Razin effect and approximate the cut-off val- 
ues po and p c by and oo, the emissivity can be calculated 
analytically (e.g. Rybicki & Lightman 1979) from Eq. (6). The 
expression for the emissivity then becomes a power law: j v oc 
y -(«-i)/2 Introducing this into Eq. (7), an analytic expression 
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which can then be compared to the results of the numerical 
code. 

In the code, the momentum integral of the emissivity, given 
in Eq. (6), is computed using the trapezoidal rule and the step 
size is refined until it reaches a desired fractional precision of 
1CT 3 . This emissivity has to be calculated at different positions 
in the wind. Then, the spatial integral of the flux, Eq. (7), is 
treated in the same way as the momentum integral, but up to 
a precision of order 1%. The number distribution is given by 
Eq. (8), while the mean synchrotron power P v is calculated us- 
ing a linear interpolation between tabulated values. This way, 
the mean synchrotron power P v needs to be calculated only 
once (as a function v/f 3 (v,p)v s (p,r)) which reduces the cal- 
culation time significantly. The next step is to approximate the 
momentum cut-off values, po and p c by and oo. We start with 
po = 1 MeV/c and p c = 10 GeV/c and then we expand the 
interval by doubling the value for p c and halving the value of 
po until the solution for the integral has converged. We can 
then compare the numerical results to the analytic solution. We 
calculated 5000 models with B t . in the range 10-100 Gauss, 
n between 1.5-7.5, 6 between 0-5 and /? max between 100- 
10000 R*. The relative errors of the numerical models com- 
pared to the analytic models did not exceed the 1 % level. 



3. Cooling mechanisms 

3.1 . Practical momentum limits 

In this section we determine practical integration limits for the 
momentum distribution, which we need to calculate the syn- 
chrotron emissivity, given by Eq. (6). 

Relativistic electrons can lose their energy through a va- 
riety of cooling mechanisms. Coulomb collisions, inverse 
Compton scattering, adiabatic cooling, synchrotron cooling 
and Bremsstrahlung cooling can all play a role. The effect 
of these cooling mechanisms has been discussed in detail by 
Chen (1992). At the high energy end of the spectrum, elec- 
trons are mainly cooled by inverse Compton scattering. This 
sets an upper limit p c to the momentum range. At the low (non- 
relativistic) end of the spectrum, Coulomb collisions are the 
dominant cooling mechanism. This sets a lower boundary po 
to the momentum range. As momenta below po or above p c 
cannot occur, the momentum range is limited to [po,p c ]- Only 
a fraction of this range actually contributes to observable radio 
emission, as can be seen from Eq. (4). For B t = 100 Gauss and 
assuming that the synchrotron radiation is formed somewhere 
between 30 and 1000/?* (see Sect. 4.2.1), the range of momenta 
that contribute to the observable radio flux is 30 - 500 MeV/c. 
It is therefore convenient to introduce an ad-hoc integration in- 
terval [pi,p2] that is contained in [po,p c ] but is large enough 
to include all momenta that produce observable radio emission. 
Within these constraints we choose [p\,p2] as small as possi- 
ble, for reasons of computational efficiency. The exact choice 
of p\ and p2 is unimportant. In the following section we shall 
derive values of p\ and p2 that can be used in the models. 



3.2. Low momentum cut-off 

Synchrotron emission at wavelengths shorter than 20 cm 
is produced by relativistic electrons with momenta above 
« 30 MeV/c. Any value /?i « 30 MeV/c can therefore be 
used as a practical lower limit to the momentum range. A con- 
venient choice of p\ is p\ = 1 MeV/c. Although the choice of 
pi does not influence the flux, it does influence the number of 
electrons (as seen in Eq. (9)). 

We still need to check whether electrons of momentum p\ 
exist, i.e. whether p\ > po (where po is the momentum cut-off 
set by Coulomb cooling). We therefore compare the acceler- 
ation time of electrons with momentum p\ = 1 MeV/c with 
the time-scale for Coulomb collisions of electrons with thermal 
ions. Collisions happen on time-scales given by (Spitzer 1956) 



to 



2 3 

mtv 3 



%ne A n % In A ' 



(12) 



where v is the speed of the relativistic electrons, n g the number 
density of the ions and In A the Coulomb logarithm. In hot-star 
winds the Coulomb logarithm is ~ 10. The time to accelerate 
an electron up to a momentum p is given by (Chen 1992) 

3 pc 

tdCC — t . / \ ~~ n> (13) 

2Au(r) eB 

where Au(r) is the shock velocity difference at radius r. For 
an electron of 1 MeV/c at a radial distance of 100 R t where the 
number density of the ions is n g = 10 7 crrT 3 and B = 0. 1 Gauss, 
we find ?d ~ 10 5 s, while f acc « 5 x 10~ 3 s (where we assume 
Au = 100 kms -1 ). Further out in the wind the ratio to/face 
increases. This means that electrons with momentum 1 MeV/c 
will not be Coulomb-cooled significantly. 

3.3. High momentum cut-off 

Synchrotron emission at wavelengths greater than 2 cm 
is produced by relativistic electrons with momenta below 
* 500 MeV/c (for B, = 100 Gauss). Any value p 2 » 
500 MeV/c can therefore be used as a practical upper limit to 
the momentum range, providing p2 < p c (where p c is the high 
momentum cut-off). 

The highest momentum attainable for relativistic electrons 
is determined by the balance between acceleration and inverse 
Compton scattering at the shock. This results in a maximum 
attainable momentum p c , that can be expressed as (Chen 1992) 



, eB(r) 
p c (r) — m e c\4n — —Au{r)r 



1/2 



, , (14) 

where cr x is the Thomson cross section and L„ the stellar bolo- 
metric luminosity. Using typical hot star values (L„ = 10 6 L©, 
v m = 2500 kms" 1 , v rot = 200 kms" 1 , R t . = 20 R and 
Au(r) = 100 km s _1 ), we can estimate a numerical value for 
the high momentum cut-off, 



p c (r) * 0.97 



GeV 



1/2 



(15) 



For Z?„ = 100 Gauss and r = 10 R», Eq. (15) gives p c = 
30 GeV/c. At larger distances from the star p c is larger. A prac- 
tical choice for p2 that satisfies p2 < p c is P2 = 15 GeV/c. 
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Table 1. Relevant stellar parameters for Cyg OB2 No. 9 adopted in 
this paper. The numbers are taken from Bieging et al. (1989, BAC), 
Leitherer (1988, L) and Herrero et al. (1999, HCVM). 



Parameter 




Ref. 


v rot sin; (km s -1 )" 


145 


BAC 


T eB (K)» 


44 500 


HCVM 


R, (R©) 


22 


HCVM 


d (kpc) 


1.82 


BAC 


v„ (km s~>) 


2 900 


L 


M (M© yr 1 ) 


2.0 x 10~ 5 


L 



" We follow White (1985) in assuming that v ro , = 250 km s~\ which 
is not in contradiction with the value for v rot sin i, 
b The wind temperature is assumed to be 0.3 times the stellar effective 
temperature (Drew 1989). 

Table 2. With the stellar parameters for Cyg OB2 No. 9 given in 
Table 1, we have the following values for the characteristic radius 
R v (see Eq. (A.7)) and for the free-free emission from the wind 
(Wright & Barlow 1975). For simplicity we assume y = 1,Z 2 = 1 
and /j. = 1.4. As a reference we also give the values of the non-thermal 
component to the flux and the observed flux with error bars. 



A 


V 




Ft 


ft 


pobs 


(cm) 


(GHz) 


(*.) 


(mJy) 


(mJy) 


(mJy) 


2 


15 


124 


1.1 


4.6 


5.7 ±0.1 


6 


5 


266 


0.6 


6.8 


7.4 ±0.1 


20 


1.4 


614 


0.3 


4.6 


4.9 ±0.1 



Now that a choice has been made for pi and p2, the inde- 
pendent parameters are reduced to B*, n, S, R mdx and /,. 

3.4. Emitting region 

We have estimated practical upper and lower limits to the mo- 
mentum range that are compatible with the high and low mo- 
mentum cut-offs due to Compton and Coulomb cooling, re- 
spectively. Of course, cooling can also change the distribution 
between these limits. Specifically, as the relativistic electrons 
move away from the shock in which they were accelerated, in- 
verse Compton scattering will rapidly cool the most energetic 
electrons. Inverse Compton scattering does not change the total 
number of relativistic particles, but quickly reduces their en- 
ergy below the energy of observable radio emission. Thus, the 
synchrotron emitting gas is concentrated in more or less nar- 
row regions behind the shocks. For the sake of simplicity we 
neglect this layered structure of the synchrotron emitting re- 
gion and assume a continuous volume of relativistic particles. 
The adequacy of this assumption is discussed in Sect. 5. 

4. Application and results 

4.1. Cyg OB2 No. 9: 21 Dec 1984 data 

We can now apply the model to Cyg OB2 No. 9 (05 If). It was 
first discovered by Abbott et al. (1984) that the spectral index 
of its radio emission does not follow the a = 0.6 law expected 



for a thermal wind source (see Fig. 1). Also, the radio emis- 
sion of Cyg OB2 No. 9 is highly variable. On one occasion 
the radio spectrum even showed a typical free-free spectrum 
(with a = 0.6). If we determine the mass loss rate from this 
radio observation, we find the same value as derived from Ha 
(Leitherer 1988). So we can assume that at that time we saw 
the underlying free-free emission from the stellar wind. Of the 
observational radio data that are available for Cyg OB2 No. 9, 
we will use the observation that has detections of emission in 
three radio wavelength bands. We chose the 21 Dec 1984 ob- 
servation with the VLA (Bieging et al. 1989) because of its 
small error bars on the detections. By fitting these observations 
we can derive the solution space of the model parameters. The 
stellar parameters adopted in our model are listed in Table 1 . 

4.2. Results 

We recall that the parameter space is spanned by n, 5, /*, 
Rmsx and B*. The parameters have a straightforward interpre- 
tation: n corresponds to a typical compression ratio in the syn- 
chrotron formation region, 6 describes the spatial dependence 
of the number distribution, /„ describes the fraction of relativis- 
tic electrons at the surface and R mdx is the outer boundary of the 
synchrotron emission region. We take B* = 100 Gauss in the 
fitting procedure, which corresponds to a local magnetic field 
offi= 1.7 xl0~ 2 Gauss atr = 500/?«. 

We calculate the synchrotron flux in the three wavelengths 
for a grid in n, 6, R mdx and /, and select those combinations of 
the parameters that fit the VLA-observations within the error 
bars. This procedure produces a strong constraint on n, 6 and 
Rm-ix, but not on /». The parameter /* is well constrained lo- 
cally (i.e. for each combination of n, 6 and R mdx ), but not glob- 
ally (i.e. for the entire range of n, 6 and R mix ). This parameter 
does play a minor role in limiting the other parameters as the 
number of relativistic electrons cannot exceed the total number 
of electrons. In principle this limitation depends on p\, but in 
practice the effect is minimal. All possible combinations of n, 
6 and R m . dx lie within the boomerang-shaped region in Fig. 2. 



The most important result from Fig. 2 is that R mdx is well con- 
strained and almost not influenced by n and 6. The outer bound- 
ary of the synchrotron emitting region must lie between 520 R t 
and 750 R*. To understand this result, it is important to realise 
that the different wavelengths play a different role in the fitting 
procedure. The fluxes at 2 and 6 cm constrain n, 6 and The 
parameters n and 6 determine the slope between 2 and 6 cm and 
/, sets the global flux level. Note that different combinations of 
n and 6 can give the same slope. The effect of R mdx on the slope 
is negligible. The flux at 20 cm is fitted by the only remain- 
ing parameter R mdx . The strong dependence of the 20 cm-flux 
on R mdx is due to the larger free-free opacity at larger wave- 
lengths. This means that more of the synchrotron emitting re- 
gion is shielded at 20 cm than at 2 and 6 cm. This also explains 
why R mdx has only a small influence on the 2 and 6 cm-fluxes. 
Note that the range in R mdx (520 - 750 /?,) extends below the 
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Fig. 2. All the combinations of n, 6 and ^ raax with B, = 100 Gauss 
that fit the VLA-observations (21 Dec 1984) for Cyg OB2 No. 9 lie 
within the boomerang-shaped region. Projections on different planes 
are plotted to situate the solutions in the parameter space. 

characteristic radio radius at 20 cm (R v (20 cm) = 614 R t ). This 
is not a contradiction as radio emission from below R v can be 
detected (see Sect. 2.3). 

While R max is insensitive to n and 6, it depends strongly 
on the magnetic field. Assuming a smaller surface magnetic 
field of 10 Gauss rather than the 100 Gauss used previously, 
we find a larger value of R mdx (see Fig. 3). R max is still well 
constrained, but now the range extends from 1 500 R t up to 
2 200 R t . The larger value of R m!lx is caused by the greater 
importance of the Razin effect at lower magnetic fields. The 
20 cm-flux is already greatly reduced by the Razin effect. To 
replenish the photons taken away by the Razin effect, the syn- 
chrotron emitting region must be larger than the value of R max 
found for a magnetic field of 100 Gauss. It is not possible to fit 
the observations with a surface magnetic field of 5 Gauss, be- 
cause essentially all radio emission is suppressed by the Razin 
effect (Sect. 2.4). In this sense the value of B, = 10 Gauss can 
be considered a lower limit for the surface magnetic field, and 
likewise the value found for R max is an upper limit to the outer 
boundary of the synchrotron emission region. 

4.2.2. n and 8 

Fig. 2 shows solutions with n extending from 2.5 to 5. Using 
the relation^ = (n + 2)/(« - 1), this translates to compres- 
sion ratios from x — 3 (moderately strong shocks) to x — 
1.75 (weak shocks). Note that strong shocks (x — 4) can- 
not explain the observations. This is consistent with what was 
found by Rauw et al. (2002), who applied a similar model to 
9 Sgr. Since time-dependent hydrodynamical simulations (e.g. 
Runacres & Owocki 2002) show a variety of shocks (and com- 
pression ratios) in the stellar wind, \ is some kind of compres- 
sion ratio average. Stronger shocks however produce more ra- 
dio emission than weaker shocks, because they accelerate more 
of the electrons into the momentum range where synchrotron 




Fig. 3. Similar figure as Fig. 2, but now calculated for a surface mag- 
netic field B, = 10 Gauss (Note the different scales compared to 
Fig. 2). 

radiation is emitted at radio wavelengths. Thus the synchrotron 
spectrum is dominated by the stronger shocks. Thus the value 
of x we obtain should be interpreted as a maximum compres- 
sion ratio in the wind, rather than a mean compression ratio. 

Note that most solutions have (5-values lower than 2. This 
means that the fraction of relativistic electrons increases further 
in the wind. This is surprising as one would expect a constant 
efficiency of the acceleration mechanism or even a decreasing 
efficiency (due to the decay of the shocks in the wind). The 
increasing relativistic fraction points to an increasing accelera- 
tion efficiency. A possible explanation is multiple acceleration 
of the electrons by subsequent shocks. Thermal electrons are 
injected in the shock front and accelerated, but electrons, al- 
ready accelerated by a previous shock, are also re-accelerated. 
This would increase the number of relativistic electrons toward 
the outer region of the stellar wind. 

Another possible explanation is that, as hydrodynamical 
simulations show a variety of shocks, some stronger shocks 
might exist in the outer parts of the wind. This would increase 
the synchrotron emission which would show up as a S < 2 in 
the parametrisation of the current model. 

4.2.3. Comparison with the Chen & White model 

We can compare our results with the results found in 
Chen & White (1991, 1994). We looked at all the models 
that fit the observations of Cyg OB2 No. 9 and derived the 
permitted range of the model parameters. We used a power 
law to describe the number distribution of relativistic elec- 
trons. Chen & White calculated the number distribution from a 
model of shocks in the stellar wind. The parameters of their 
synchrotron model which reproduce the radio spectrum of 
Cyg OB2 No. 9 were n = 2.6, 6 < 2 and fl max = 500fl«. Even 
though we used different methods to describe the momentum 
distribution of relativistic electrons in the wind, the parameters 
of their synchrotron model are very similar and can be found 
near the solution space of our model. 
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5. Discussion and conclusions 

In this paper we have presented a theoretical model for 
non-thermal radio emission from single hot stars, in terms 
of synchrotron emission from electrons accelerated in wind- 
embedded shocks. Applied to a specific observation of 
Cyg OB2 No. 9, the model produces constraints on almost all 
of its parameters (/„ being the only exception). A major result 
is that the outer boundary of the synchrotron emitting region 
(^max) is well constrained (520 - 750 R* for B„ = 100 Gauss). 
As the efficiency of Compton cooling prevents relativistic elec- 
trons from carrying their kinetic energy far from the shock, this 
means that our models predict that shocks persist to at least 
500 R*. This is consistent with recent hydrodynamical calcula- 
tions (Runacres & Owocki 2003). 

Another interesting result is the absence of very strong 
shocks (those with compression ratio \ near 4) in the outer 
wind. We find shocks with x - 3 at most. This is consis- 
tent with the hydrodynamical prediction that shocks become 
weaker as they move out with the stellar wind. 

Within the range of solutions we find both 6 ~ 2 and 6 < 2. 
The physical picture associated with these different values of 6 
is rather different too. The S » 2 solutions point to a constant 
fraction of relativistic electrons. It can be seen from Fig. 2 that 
they also correspond to weak shocks (x ^ 2). This offers a 
natural explanation for R max : the shocks weaken as they move 
out with the flow and beyond R m . dx they are to weak to produce 
significant synchrotron emission. 

However, most of the solutions have 6 < 2, which points 
to the relativistic electron fraction increasing outwards. This 
could indicate an increasing compression ratio of the shocks, 
or at least the presence of a stronger shock in the outer part 
of the synchrotron emitting region. While this could inciden- 
tally have been the case at the time of observation, it does not 
sit well with the hydrodynamical prediction that shocks decay 
as they move outward. It is also hard to imagine why these 
shocks would rather suddenly disappear when they reach 7? max . 
A more plausible explanation of 6 < 2 is multiple accelera- 
tion: particles in the outer wind have had greater opportunity to 
pass through multiple shocks than particles in the inner wind. 
The 5 < 2 solutions are also those with larger compression 
ratios. In the case of multiple acceleration the momentum dis- 
tribution is no longer a pure power law ~ p~ n , but a flatter 
function (White 1985). Simply put, this means that the mo- 
mentum distribution for multiple acceleration by a number of 
weaker shocks resembles a distribution corresponding to a sin- 
gle stronger shock. This is then again consistent with the idea 
of an outer boundary of the synchrotron emitting region, as all 
the shocks involved in the acceleration are weak. 

In this paper we have assumed a power-law momentum dis- 
tribution with a single compression ratio, where all of the vol- 
ume can contribute to the emission. The efficiency of inverse- 
Compton cooling, however, will limit the synchrotron emission 
region to more or less narrow layers behind the shocks. The pri- 
mary effect of this is to reduce the number of electrons in the 
momentum range relevant for radio emission. This means that 
the values of f, we find in our fitting procedure are underesti- 
mated as it is this parameter that sets the overall level of syn- 



chrotron emission, and must mimic the quenching of electrons 
by inverse Compton scattering. Because /» is not well con- 
strained by the adopted fitting procedure, this is of no conse- 
quence for our results. As a secondary effect, inverse-Compton 
cooling can also change the momentum distribution, as it af- 
fects the high momentum electrons more than the low momen- 
tum electrons. The importance of this effect will be investigated 
in a subsequent paper. 

Hydrodynamical models predict a variety of compression 
ratios. Due to the strong dependence of the emission on the 
compression ratio, the strongest shock will dominate. In that 
case, the radiation would not come from the whole volume, but 
only from behind the strongest shocks. This may also explain 
the variability that is a characteristic feature of non-thermal ra- 
dio emission (Bieging et al. 1989). As the stronger shocks move 
in and out of the synchrotron emitting region they lead to vari- 
ability. 

The presence of relativistic electrons in the wind will also 
lead to non-thermal X-rays, due to the inverse Compton process 
(Chen & White 1991). The non-thermal X-rays should show 
up as a hard X-ray tail in a spectrum. Clear evidence for such 
a non-thermal tail has remained elusive. Rauw et al. (2002) 
show that the hard X-ray tail of the non-thermal radio emit- 
ter 9 Sgr does not have the expected power-law dependence 
for non-thermal X-ray emission. The tail can also be modelled 
with a high-temperature (> 2x 10 7 K) thermal plasma. It should 
be noted that the non-thermal X-ray emission is formed much 
closer to the stellar surface (where there is an abundant sup- 
ply of UV photons) than the non-thermal radio emission. One 
could therefore surmise that the shocks closer to the star are 
weaker and therefore result in an insufficient number of rela- 
tivistic electrons to create detectable non-thermal X-rays. This 
suggestion, however, does not appear to be consistent with cur- 
rent hydrodynamical models, as the shocks are strongest in the 
inner wind and decay as they move outward. The hard X-ray 
tail of 9 Sgr could also be explained by colliding wind emis- 
sion from a long-period binary companion (Rauw et al. 2002). 

To obtain non-thermal radio emission, all that is required 
is a magnetic field and shocks in the wind. As both should be 
present in any early-type star, the question arises why not all 
these stars are non-thermal emitters. One possible answer is 
that the magnetic field in thermal stars is lower. A lower mag- 
netic field means that less synchrotron radiation is generated 
and that the Razin effect is more efficient at suppressing it. This 
explanation suggests that searches to detect magnetic fields in 
O-type stars might have a better chance in non-thermal radio 
emitters. Another possible explanation is that the synchrotron 
emitting region (essentially determined by the characteristic 
free-free radius and R max ) is too small, so that the synchrotron 
radiation emitted falls below detection levels. 

Finally, one needs to ask the uncomfortable question 
whether there is any such thing as synchrotron emission 
from single O stars. For Wolf-Rayet stars, the connection be- 
tween non-thermal emission and binarity is firmly established 
(Dougherty & Williams 2000). For O stars, the situation is less 
clear. Roughly half of the non-thermal O stars are not known to 
be members of a binary system. However, the truism that ab- 
sence of evidence is not evidence of absence obviously holds. 
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Also, there is an indication that one of the archetypal sin- 
gle non-thermal O star, 9 Sgr, might actually be in a binary 
(Rauw et al. 2002). The possibility that all non-thermal emit- 
ting O stars are binaries cannot be excluded. 
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Appendix A: Flux: exact and approximated integral 

We calculate the flux following Wright & Barlow (1975) for 
the free-free emission, but we replace xB v , where \ is the free- 
free absorption coefficient and B v the source function (e.g. the 
Planck function), by the synchrotron emissivity j v . The exact 
integral for the flux is 

F v = — d^j dlj v [r(l)K TAq ''\ (A.l) 

where D is the distance to the star, q the impact parameter, / the 
distance along the line of sight with an observer at I = -co and 
R m . dx the outer boundary for the synchrotron emission region. 
The optical depth T v (q, I) is defined as 



r v {q, I) 



/ 

■J -ex 



61' 



K s yA 2 
(q 2 + l' 2 ) 2 ' 



(A.2) 



with A = M /(4n //nigv m ), where the symbols have their usual 
meaning, and Ks the free-free absorption coefficient given by 
Allen (1973). Note that by using as the lower integral bound- 
ary in Eq. (A.l), we in fact replaced the star by stellar wind 
material. This mathematical simplification has no consequence 
as the free-free optical depth is so large that any contribution of 
the star can be neglected anyway. 

The equation for r v (q,l) can be integrated to (using 
Abramowitz & Stegun 1965, 3.3.24) 



r Y (q, I) = 



K B yA 2 
2q 3 



l/q 



1 + 



+ arctan 



qj 2 



(A.3) 



Introducing polar coordinates, 

q = rsinO 
I = rcos6 
dqdl = rdrdO 
we find 

Kay A 2 

T v (q,l) = — ^-r- (sinfl cos 9- 6 + n) 
2r 3 sin 8 

and for the flux 
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with R v the characteristic radius 
(Wright & Barlow 1975) expressed as 

4 



rQ)(i) 
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(A.5) 

(A.6) 
of emission 

(A.7) 
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Fig. A.l. The full line is the geometric function G(r/R v ), while the 
dashed line represents the Heaviside function. 



and 

G(. 



d#sin6> 



(A.8) 



x exp 



(r(3»V 

x 3 sin 3 6 



(sinflcosfl - 6 + n) 



The 6 integration is only from to n, because the original in- 
tegral only considered positive q values (a factor of 2 was in- 
cluded in the original integration to compensate for this). The 
function G(r/R v ) is plotted on Fig. A.l. Note that G(r/R v ) can 
be approximated by a Heaviside function: H = if r < R Y 
and H — 1 if r > R v . This brings us back to the approxima- 
tion Eq. (7) for the flux. The function G(r/R y ) needs to be cal- 
culated only once, so that evaluating Eq. (A.6) does not take 
significantly more time than Eq. (7). 
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